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Abstract: Climate change may cause shifts in the natural range of species especially for those that are 
geographically restricted and/or endemic species. In this study, the spatial distribution of five endemic 
and threatened species belonging to the genus Onosma (including O. asperrima, O. bisotunensis, O. kotschyi, O. 
platyphylla, and O. straussii) was investigated under present and future climate change scenarios: RCP2.6 
(RCP, representative concentration pathway; optimistic scenario) and RCP8.5 (pessimistic scenario) for the 
years 2050 and 2080 in Iran. Analysis was conducted using the maximum entropy (MaxEnt) model to 
provide a basis for the protection and conservation of these species. Seven environmental variables 
including aspect, depth of soil, silt content, slope, annual precipitation, minimum temperature of the 
coldest month, and annual temperature range were used as main predictors in this study. The model 
output for the potential habitat suitability of the studied species showed acceptable performance for all 
species (i.e., the area under the curve (AUC)>0.800). According to the models generated by MaxEnt, the 
potential current patterns of the species were consistent with the observed areas of distributions. The 
projected climate maps under optimistic and pessimistic scenarios (RCP2.6 and RCP8.5, respectively) of 
2050 and 2080 resulted in reductions and expansions as well as positive range changes for all species in 
comparison to their current predicted distributions. Among all species, O. bésotunensis showed the most 
significant and highest increase under the pessimistic scenario of 2050 and 2080. Finally, the results of 
this study revealed that the studied plant species have shown an acute adaptability to environmental 
changes. The results can provide useful information to managers to apply appropriate strategies for the 
management and conservation of these valuable Iranian medicinal and threatened plant species in the 
future. 
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1 Introduction 


Anthropogenic activities have increased greenhouse gases in the atmosphere and have been 
causing climate change most especially since the mid-20" century (Sala et al., 2000; Pereira et al., 
2010; Peñuelas et al., 2013). Global warming affects all levels of biodiversity from individuals to 
ecosystems (Parmesan, 2006; Bellard et al., 2012). The most important effects are increases of 
global temperatures, variations in the lifecycles of plants and animals, changes in the species 
distribution, loss of biodiversity, drought, wildfires, and others (Amedie, 2013). The geographic 
distribution pattern of species depends on many environmental variables and the ability of 
individual species to adapt to the new conditions. Climate change is one of the most important 
factors known to alter species distributions by reducing or expanding their geographical range 
(Tallis et al., 2008; Garcia et al., 2014; Alamgir et al., 2015). The Intergovernmental Panel on 
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Climate Change (IPCC) estimates a high risk of extinction for plant species in the 21* century 
and further projects that about 58% of plants will lose their ecological niche and natural habitats 
to climate change in 2080 (Warren et al., 2013). Understanding the current and future spatial 
response patterns of plant species to climate change is a basic issue in ecology and conservation 
biogeography (Guisan et al., 2013; Guillera-Arroita et al., 2015). 

Iran with more than 8100 vascular plant species has a diverse flora among the countries of 
southwestern Asia (Noroozi et al., 2016). A wide array of climatic conditions and pedological 
diversity have led to a high level of species richness in Iran (Hedge and Wendelbo, 1978) and 
approximately 24%—30% of Iranian vascular plants species are endemics (Akhani, 2006; Noroozi 
et al., 2016). Over the past decades, however, land use changes, overgrazing, overharvesting, as 
well as invasive and newly-introduced species have been threatening the plant diversity 
(Mehrabian, 2015; Sayadi and Mehrabian, 2016; Noroozi et al., 2018; Mehrabian et al., 20202). 
Further, it is expected that climate change will have an additional impact on this extensive 
diversity (Parmesan, 2006; Loarie et al., 2008; Abdelaal et al., 2019) as it is estimated that if the 
CO» concentration doubles by the year 2100, the average temperature in Iran will increase by 
1.5°C—4.5°C (Roshan et al., 2011). 

As endemic species are restricted to a specific geographic area and have unique genetic 
reserves (Bonn et al., 2002; Fois et al., 2018; Abdelaal et al., 2019), climate change may cause 
shifts in the natural range of those species unable to adapt to new climatic conditions and 
consequently may face a higher risk of endangerment or even extinction (Loarie et al., 2008; 
Abdelaal et al., 2019; Bender et al., 2019). The genus Onosma is one of the most important 
genera of Boraginaceae in Iran as it is comprised of a large number of species and has a high rate 
of endemism. So far, 54 species of this genus have been reported in Iran, 24 of which are endemic 
to the flora of Iran (Khatamsaz, 2002; Mehrabian and Amini, 2018). Moreover, this genus has 
medicinal properties and the root of its plants contain a major active chemical component known 
as shikonin (Sut et al., 2017). Shikonin derivatives have antibacterial, anti-inflammatory, and 
anti-tumor properties, and could bypass cancer drug resistance (He, 2009; Liu et al., 2010; Noula 
et al., 2010). 

The limited distribution of endemic and rare species of this valuable genus as well as 
destructive human activities and droughts in recent years, have caused such severe damage to so 
many species that about 50% are now found on the International Union of Conservation of 
Nature's (IUCN) Red List of Threatened Species (Mehrabian, 2015). The mounting destruction 
motivates us to investigate the repercussions of climate change on the future spatial distributions 
of endemic and threatened species of this genus, including O. asperrima Bornm (NT, near 
threatened), O. bisotunensis Attar (EN, endangered), O. kotschyi Boiss. (NT), O. platyphylla H. 
Riedl (NT), and O. straussii (Riedl) Khat. (NT) (Fig. 1). Climate change may negatively intensify 
the severity of the impact on the survival of these species. It is also critical to define the areas that 
these species with narrow niches inhabit or most likely inhabit in order to apply an appropriate 
strategy for their conservation (Dubuis et al., 2011; Kaky and Gilbert, 2016). 

In this regard, species distribution modeling (SDM) has become a key method in ecology and 
conservation biogeography to predict the distribution of a species across geographic space and 
time using environmental data (Margules and Pressey, 2000; Groves et al., 2002; Graham et al., 
2004; Peterson and Soberón, 2012). Numerous algorithms using presence and/or absence data 
have been developed to predict the geographical distribution of a given species (Soberón and 
Peterson, 2005; Elith et al., 2006; Elith and Leathwick, 2007, 2009). The maximum entropy 
(MaxEnt) modeling, however, is a machine learning algorithm that has been used extensively for 
distribution modeling since it is recognized as one of the best performing methods for modeling 
species presence-only data (Elith et al., 2006). MaxEnt is also recognized as a high-performance 
algorithm for predicting species distributions with few occurrence presence points as endemic 
species generally have low distributions (Elith et al., 2006; Phillips et al., 2006). Many ecological 
studies have used the SDM approach throughout the world (Ródder and Weinsheimer, 2009; 
Aragón et al., 2010; Rubidge et al., 2011; Khanum et al., 2013; Kujala et al., 2013; Legault et al., 
2013; Adams-Hosking et al., 2015; Bleyhl et al., 2015; Luo et al., 2015; Sen et al., 2016; Ulrey et 
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al., 2016). Despite the very rich species diversity in Iran, however, there have been very few 
studies in this field regarding local flora (e.g., Ardestani et al., 2015; Mazangi et al., 2016; 
Abolmaali et al., 2018). Hence, the objective of our study is to predict the current and future 
spatial distributions of O. asperrima, O. bisotunensis, O. kotschyi, O. platyphylla, and O. straussii 
for the first time to forecast the enlargement and direction of niche shifts in scenarios projecting 
climate change (RCP2.6 (RCP, representative concentration pathway) and RCP8.5) for the years 
2050 and 2080 using a set of environmental variables and the MaxEnt model. 


Fig. 1 Photographs of O. Asperrima (a), O. Platyphylla (b), O. Straussii (c), O. Kotschyi (d), and O. 
Bisotunensis (e) in Iran 


2 Materials and methods 


2.1 Study area 


The study was conducted in Iran, and this country consists of the vast Iranian Plateau and has an 
area of about 1.65x10° km?. As the second largest country in the Middle East, Iran is located in 
the dry belt of Asia. It shares borders with Armenia and the Republic of Azerbaijan in the 
northwest, the Caspian Sea in the north, Turkmenistan in the northeast, Afghanistan and Pakistan 
in the east, Turkey and Iraq in the west, and the Persian Gulf and the Gulf of Oman in the south. 
Climatic factors heavily influence the distribution patterns of plant taxa, while other ecological 
factors have more localized effects (Zahran, 2010). Orogenic massifs surround the central plateau 
of Iran and act as a natural barrier to the penetration of humidity to these zones. This also shapes 
the diverse rainfall patterns. Iran is one of the most mountainous countries in the world. Its 
landscape is surrounded by several high, rugged mountain ranges such as the Caucasus, Zagros, 
Alborz, and Kopet-Dagh. The Zagros system has a northwest-southeast orientation extending 
from eastern Turkey to southwestern Iran (Fisher, 1968). It forms a wall among the Iranian 
Plateau, the Mesopotamian, and the Persian Gulf. It has a mean elevation of 1200 m a.s.l. and the 
highest peak, Zard Kuh, reaches 4231 m a.s.l. (Homke et al., 2004). The study covers the 
Mediterranean macro-bioclimate (Rivas-Martínez et al., 1997, 1999). 


2.2 Species occurrence data 


We have selected five species (O. asperrima, O. bisotunensis, O. kotschyi, O. platyphylla, and O. 
straussii) of Onosma genus because enough data were available for modeling. Distributional data 
on characterizing species were gathered from (i) field surveys during 2009—2017; (ii) literature 
records available in Flora Iranica (Riedl, 1967) and Flora of Iran (Khatamsaz, 2002); and (iii) 
historical data available in the herbaria of HSBU, IRAN, TARI, WU, BASU, and B (the 
herbarium acronyms follow the study of Thiers (2019)). Due to the lack of careful and reliable 
absence species distribution data, only presence data were used in this study. Distribution map 
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(presence data) of O. asperrima, O. bisotunensis, O. kotschyi, O. platyphylla, and O. straussii in 
Iran is shown in Figure 2. The point inputs to the models developed in this study were collected 
from their habitats in the west and southwest of Iran such as Fars, Lorestan, Kurdistan, 
Khuzestan, Kermanshah, Hamedan, Markazi, Ilam, and Kohgiluyeh and Boyer-Ahmad provinces 
and certain areas in the south and central provinces like Kerman and Isfahan. 
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Fig. 2 Distribution map (presence data) of O. asperrima, O. bisotunensis, O. kotschyi, O. platyphylla, and O. 
straussii in Iran 


2.3 Selection of environmental variables 


In this study, according to expert opinion and the ecological need of the species, we initially used a 
total of 13 environmental variables related to the distribution of Onosma as follows: BIO1 (annual 
mean temperature), BIOS (maximum temperature of warmest month), BIO6 (minimum 
temperature of coldest month), BIO7 (annual temperature range (BIOS—BIO6)), BIO12 (annual 
precipitation), aspect, elevation, slope, solar radiation per month, land use, silt content, soil organic 
carbon content, and depth of soil. Afterwards, collinearity among environmental variables was 
tested by 
one of them was excluded in order to avoid co-linearity (Elith et al., 2010). The source of selected 
variables after correlation test is available in Table 1. To represent climate change influences, we 
used projected future climate variables for 2050 and 2080 with empirically downscaled bioclimatic 
data downloaded from the CCAFS website (Climate Change, Agriculture and Food Security; 
http://www.ccafs-climate.org), and the average of 16 General Circulation Models (GCMs) under 
optimistic (RCP2.6) and pessimistic (RCP8.5) greenhouse-gas emissions scenarios. The resolution 
of environmental variables used in the study was 30 arc-seconds (ca. 1 kmx1 km). 


2.4 Modeling process and evaluation 


MaxEnt model (Phillips et al., 2006) was utilized for modeling the current and future habitat 
suitability of species. MaxEnt (jar file v3.4.1) was utilized through the dismo package v1.1-4 
(https://rspatial.org/raster/sdm/) in the R v3.2.3 programming environment (R Core Team, 2018). 
According to some previous studies (Bosso et al., 2013; Vasconcelos et al., 2014; Fois et al., 
2018), we predominantly used the MaxEnt model when the data points included presence-only 


Table 1 Source of selected variables after correlation test and estimates of their permutation importance in 
maximum entropy (MaxEnt) modeling for the studied Onosma species 


Permutation importance (46) 


Source Variable O. O. O. O. O. 
asperrima — bisotunensis kotschyi platyphylla — straussii 


Bioclimatic variables BIO6 1.8 0.1 1.6 5.4 4.5 
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(www.worldclim.org) BIO7 0.3 37.7 7.0 24.2 32.5 

BIO12 11.9 34.9 7.6 ila 28.4 

Topographic variables Aspect 0.7 2.7 0.5 17 0.7 
(www.worldgrids.org) Slope 70.7 DS 754 35.7 12.7 
Edaphic variables (www.soilgrid.org; Silt content 42 0.9 7.5 23.1 12.0 
www.isric.org) Depth of soil 10.4 21.3 0.5 8.2 9.3 


Note: BIO6, minimum temperature of coldest month; BIO5, maximum temperature of warmest month; BIO7, annual temperature range; 
BIOI2, annual precipitation. BIO7=BIO5—BIO6. Aspect and slope were all derived from elevation. 


with a limited number of records. This is because the MaxEnt fits models with varying levels of 
complexity with respect to the amount of available data. Thus, for the species with few presence 
records, it fits simpler models (Elith et al., 2011). The models were evaluated using 10-fold 
cross-validation. In cross-validation, data were randomly divided into ten parts; nine parts were 
used for model fitting and the fitted model was then used to evaluate the holdout part (Valavi et 
al., 2019). For predicting the current and future suitability, a single MaxEnt model with full 
dataset was re-fitted. 

We also considered permutation importance in order to define main environmental variables 
which have influenced the potential distributions of the studied species (Abdelaal et al., 2019). To 
assess the accuracy of the modeling results, we computed the area under the curve (AUC) of the 
receiver operating characteristic curve (Lobo et al., 2008). The AUC score is a powerful tool for 
measuring model performance because of its independence from threshold selection (Yi et al., 
2016; Fois et al., 2018). The AUC shows the power of the model to discriminate presences from 
random background (Phillips et al., 2009) The AUC ranges between 0.000 and 1.000, with 0.500 
showing a random prediction performance and 1.000 indicating a perfect discrimination. Values 
under 0.500 indicate models worse than random (Elith et al., 2006). The equal training sensitivity 
and specificity of MaxEnt output was considered as proper threshold for the model prediction 
(Liu et al., 2013). 


3 Results 


3.1 | Model performance and the key environmental factors influencing species 
distributions 


After the correlation test, we selected seven environmental variables for modeling (Table 1). The 
modeling outputs for the potential habitat suitability of O. asperrima, O. bisotunensis, O. 
kotschyi, O. platyphylla, and O. straussii showed a perfect predictive performance with AUC 
values higher than 0.800 (0.823, 0.973, 0.922, 0.943, and 0.895, respectively). Considering 
permutation importance, slope, BIO12, and depth of soil were the main environmental variables 
to have influenced the potential distribution of O. asperrima (Table 1). In this regard, BIO7, 
BIO12, and depth of soil were most important for O. bisotunensis. For O. kotschyi, most 
important variables were slope, BIO12, and silt content. For O. platyphylla, variables of slope, 
BIO7, and silt content were important; and finally for O. straussii, variables of BIO7, BIO12, and 
slope were important (Table 1). As can be further seen in Table 1, BIO6 and aspect proved less 
important than other variables for all species. 


3.2 Potential habitat suitability of the studied species in current and future conditions 


Species distribution maps showed that the potential suitable habitats of O. asperrima and O. 
kotschyi are currently located in some parts of the Kurdo-Zagrosian and Fars-Kerman 
phytogeographic sub-provinces in Lorestan, Chaharmahal and Bakhtiari, Kohgiluyeh and 
Boyer-Ahmad, Fars, and Kerman provinces (Figs. 3 and 4). The MaxEnt model predicted that the 
potential suitable habitats of O. bisotunensis lie in the northern parts of the Kurdo-Zagrosian 
phytogeographic zone in Kermanshah and Kurdistan provinces as well as in some parts of West 
Azarbaijan (Fig. 5). The suitable habitats of O. platyphylla are mainly located in the central 
Zagros (Kermanshah, Kurdistan, Lorestan, Chaharmahal and Bakhtiari, and Kohgiluyeh and 


202101.00050v1 


chinaXiv 


JOURNAL OF ARID LAND 


ChinaXiv& ERAT 


Boyer-Ahmad) as well as in other parts of the Kurdo-Zagrosian phytogeographic sub-province 
(West Azarbaijan, Fars, and Kerman provinces) (Fig. 6). 
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Fig. 3 Projection of potential distribution maps of O. asperrima currently and under climate change scenarios 
RCP2.6 and RCP8.5 in 2050 and 2080. (a), current prediction; (b), future prediction under RCP2.6 in 2050; (c), 
future prediction under RCP2.6 in 2080; (d), future prediction under RCP8.5 in 2050; (e), future prediction under 
RCP8.5 in 2080. 
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Fig. 4 Projection of potential distribution maps of O. kotschyi currently and under climate change scenarios 
RCP2.6 and RCP8.5 in 2050 and 2080. (a), current prediction; (b), future prediction under RCP2.6 in 2050; (c), 
future prediction under RCP2.6 in 2080; (d), future prediction under RCP8.5 in 2050; (e), future prediction under 
RCP8.5 in 2080. 
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Fig. 5 Projection of potential distribution maps of O. bisotunensis currently and under climate change scenarios 
RCP2.6 and RCP8.5 in 2050 and 2080. (a), current prediction; (b), future prediction under RCP2.6 in 2050; (c), 
future prediction under RCP2.6 in 2080; (d), future prediction under RCP8.5 in 2050; (e), future prediction under 
RCP8.5 in 2080. 
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Fig. 6 Projection of potential distribution maps of O. platyphylla currently and under climate change scenarios 
RCP2.6 and RCP8.5 in 2050 and 2080. (a), current prediction; (b), future prediction under RCP2.6 in 2050; (c), 
future prediction under RCP2.6 in 2080; (d), future prediction under RCP8.5 in 2050; (e), future prediction under 
RCPS.5 in 2080. 


The potential suitable habitats of O. straussii are located in the Kurdo-Zagrosian 
phytogeographic sub-province and some regions of the central phytogeographic sub-provinces in 
Kermanshah, Lorestan, Markazi, Hamedan, and Khuzestan provinces (Fig. 7). The projected 
climate maps under optimistic and pessimistic scenarios (RCP2.6 and RCP8.5, respectively) of 
2050 and 2080 resulted in reductions and expansions as well as positive range changes for all 
species in comparison to their current predicted distributions (Figs. 3—7; Tables 2 and 3). Among 
all species, O. bisotunensis showed the most significant and highest increase under the 
pessimistic scenarios of 2050 and 2080. Among all, O. straussii exhibited the highest increase 
under the optimistic scenarios of 2050 and 2080. 
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Fig. 7 Projection of potential distribution maps of O. straussii currently and under climate change scenarios 
RCP2.6 and RCP8.5 in 2050 and 2080. (a), current prediction; (b), future prediction under RCP2.6 in 2050; (c), 
future prediction under RCP2.6 in 2080; (d), future prediction under RCP8.5 in 2050; (e), future prediction under 


RCPS.5 in 2080. 


Table 2 Percentages of gain, loss, and change in the distributions of the studied species under climate change 
scenario RCP2.6 in 2050 and 2080 


Climate change scenario RCP2.6 


Species 2050 2080 
Gain (96) Loss (96) Change (96) Gain (96) Loss (96) Change (96) 
O. asperrima 42.80 9.43 35.9 38.67 11.07 27.60 
O. bisotunensis 50.01 27.78 2223 54.74 27.90 26.84 
O. kotschyi 42.47 15.05 27.42 38.87 18.02 20.86 
O. platyphylla 48.57 16.22 32.30 47.71 18.38 2933 
O. straussii 64.37 14.71 49.66 65.38 14.13 51.76 


straussii. 


Note: Threshold was 0.369 for O. asperrima, 0.654 for O. bisotunensis, 0.349 for O. kotschyi, 0.398 for O. platyphylla, and 0.397 for O. 


Table 3 Percentages of gain, loss, and change in the distributions of the studied species under climate change 
scenario RCP 8.5 in 2050 and 2080 


Climate change scenario RCP8.5 


Species 2050 2080 
Gain (%) Loss (%) Change (%) Gain (%) Loss (%) Change (%) 
O. asperrima 58.83 7.50 46.33 52.03 8.91 43.12 
O. bisotunensis 13:07) 23.51 50.46 164.34 9.11 155.23 
O. kotschyi 51.16 11.78 39.38 49.01 13.03 35.98 
O. platyphylla 46.79 iaa 29.07 45.46 PPAS 23.06 
O. straussii 65.17 14.52 51.20 69.65 14.58 55.07 


straussii. 


4 Discussion 


Note: Threshold was 0.369 for O. asperrima, 0.654 for O. bisotunensis, 0.349 for O. kotschyi, 0.398 for O. platyphylla, and 0.397 for O. 


This study has demonstrated the potential geographical distributions of five endemic species of 
the Onosma in Iran under both current and future climate scenarios. The current distribution 
patterns of these taxa are often concentrated in the western parts of the country. The Zagros 
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Mountains are one of the most important centers of speciation in Iran, with a rich and diverse 
flora and fauna. These mountain ranges are located in the transitional zones between the 
Irano-Turanian and Mediterranean areas, which have created a variety of habitats and ecosystems 
for many plant species with a high degree of endemism (Zohary, 1973; Hedge and Wendelbo, 
1978; Mehrabian, 2015; Mehrabian et al., 2020b). The phyto-geographical studies of this genus 
indicate that the western slopes of Zagros Mountains are one of the important plant areas and a 
priority for conservation of Onosma species in Iran (Mehrabian, 2015). Seventy-five percent of 
all endemic species of the genus including O. asperrima, O. bisotunensis, O. kotschyi, O. 
platyphylla, and O. straussii are distributed in the central and northern Zagros mountain 
ecosystem. Zagros is one of the most important biodiversity zones in the world and one of the 
most vital endemic centers of this genus (Mehrabian, 2015). Many endemic and crucial species of 
this genus are endangered in Zagros due to their limited distributions and the wide range of 
threats they face. Unfortunately, these ecosystems have been at risk of destruction for various 
reasons. Land conversion, overgrazing, and fire have reduced the number of valuable Zagrosian 
plant species in recent decades. In addition, several studies have shown the negative effects of 
climate change on the Zagros mountain ecosystem (Hosseini and Asghari, 2012; Valavi et al., 
2018). Due to the fact that the O. asperrima, O. bisotunensis, O. kotschyi, O. platyphylla, and O. 
straussii species are endangered (EN) according to the IUCN Red List of Threatened Species 
(Mehrabian, 2015) and the suitable habitats of this species are located in the Zagros Mountains, 
their habitats must be protected accordingly. 

The results obtained in the present study showed that BIO12, BIO7, silt content, slope, and 
depth of soil are generally key to the geographic distributions of the studied species. Climate is a 
dominant component in ecology and affects the growth of plants and their distribution patterns. 
Therefore, climate variables act as a natural selection factor in the formation of diverse vegetation 
(Adams, 2007). Several ecological studies (Mehrabian, 2015; Moradi et al., 2019) in different 
bioclimatic zones of Iran emphasize the importance of climate in the establishment of endemism, 
diversity, and genetic divergence in Onosma. For example, precipitation (BIO12) is actually a key 
eco-factor for the habitat suitability of most species in the present study. By affecting seed 
germination (Quevedo-Robledo, 2010), seedling growth and survival (Padilla and Pugnaire, 
2007), and phenology (Mathias and Chesson, 2013), precipitation plays a determinant role in 
species richness, species distribution patterns, and diversification of plant species (Pausas and 
Austin, 2001; Yan et al., 2015). This environmental factor has a fundamental role in the potential 
distribution models of certain Iranian plant species including Astragalus caragana Fischer & C. 
A. Meyer (Ardestani et al., 2015) and Daphne mucronata Royle (Abolmaali et al., 2018). Further, 
the temperature variable such as BIO7 is a main environmental factor influencing the spatial 
distributions of most species in this study. 

The distributions of plant taxa also mainly depend on edaphic factors (silt content and depth of 
soil) and geomorphological factors (slope and aspect) (Hanson and Churchill, 1962; Guisan and 
Thuiller, 2005). The edaphic parameters, unlike climate variables, are not homogeneous across a 
landscape and can vary based on the type of parental materials (Anderson, 1988), 
geomorphologic situations (Ceddia et al., 2009), and land use (Mwanjalolo Jackson-Gilbert, 
2015) from place to place. Therefore, they should always be considered in SDM studies (Velazco 
et al., 2017). Several studies (Mehrabian, 2015; Sayadi et al., 2017; Moradi et al., 2019) on the 
ecology of this genus in Iran emphasize the role of soil factors such as depth and silt content in 
the diversity and endemism of this genus. Geomorphology is also a main ecological factor 
influencing the distributions and diversity of vegetation and plant taxa in mountainous zones 
(Cantlon, 1953; Coblentz and Riitters, 2004). Accordingly, it controls spatial redistribution of 
sunlight, heat, water, and soil nutrients (Parker and Branner, 1982; Feng et al., 2011). Several 
studies on plant diversity in the Taihang Mountains of China (Li and Zhang, 2006), species 
diversity of the forest communities in the southern Taihang Mountains of China (Ru et al., 2006), 
topographic factors on vegetation mosaics and tree diversity in the Chihuahuan Desert of North 
America (Poulos and Camp, 2010), and vegetation and plant species diversity in the southern 
slopes of Vitosha Mountain, Southeast Bulgaria (Dyakov and Nikolay, 2014), emphasize the 
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importance of slope in species diversity and distribution patterns in plant taxa. In addition, 
Mehrabian (2015) and Moradi et al. (2019) have confirmed these results in the context of 
ecological studies on Onosma in Iran. 

Species will have three responses to climate change: adapting to the new climate, migrating to 
a suitable habitat, and going extinct (Yousefi et al., 2020). Climate change will probably affect 
most species negatively but some may benefit if the amount of suitable habitat increases. 
Identifying winners and losers is becoming an important topic in climate change studies (Muths et 
al., 2017; Kafash et al., 2018; Yousefi et al., 2020). Having identified winners and losers it would 
be possible to manage each species according to their specific response to changes. In this study, 
we modelled the impacts of climate change on the spatial distributions of five endemic species 
under optimistic and pessimistic scenarios (RCP2.6 and RCP4.5, respectively) of 2050 and 2080. 
Our results showed that all species will respond similarly to changes. As a matter of fact, suitable 
habitat for all five species will increase in the future and these species will benefit from a change 
in climate, even though they are endemic and threatened. Our results are in line with previous 
studies which predicted range expansions (Kafash et al., 2016; Farashi and Erfani, 2018; Kafash 
et al., 2018) for different taxonomic groups in the country. 

The studied plant species have shown acute adaptability to environmental changes. On the 
basis of the physiology, genetics, and ecological features of different taxa, climate change will 
have a variety of effects on their future distribution patterns. Onosma is centered mainly in 
xero-habitats, including rock and sandy soils as well as serpentine geological formations (Cecchi 
et al., 2011) in the mountainous habitats of the Irano—Turanian regions of Asia along with the 
Mediterranean region, especially in Iran and Turkey. All of the above ecological conditions 
indicate the high adaptability of this species to stressful environmental conditions, especially to 
hot, dry habitats. This was confirmed by numerous geo-botanical studies (Mehrabian, 2015; 
Moradi et al., 2019) that have been conducted in Iran. In addition, several endemic species of 
Onosma (e.g., O. straussii, O. kotschyi, O. bistounensis, O. mozaffariani, O. sheidai, O. 
chrysocaheta, O. sarvestanica, and others) grow in the xeric steppe habitats of Iran; Mehrabian 
(2015) has proven that they are well adapted to these stressful habitats. All of the above reasons 
reinforce the possibility of habitat expansion of these species under future climate change 
scenarios. Similarly, some ecological modeling on drought-friendly species such as Capparis 
spinosa L. (Ashraf et al., 2018), Ambrosia artemisiifolia L. (Adhikari et al., 2019), Ambrosia 
trifida L. (Adhikari et al., 2019), and Solanum carolinense L. (Adhikari et al. 2019) confirmed 
these results. Additionally, climate change can alter the competitive interactions and species 
composition within communities (Howden et al., 2003). Therefore, a reduction in the competition 
of its companion species can lead to the expansion of these resistant species. On the basis of field 
observations, Mehrabian (2015) revealed that O. strausii rapidly spreads due to reduced 
competition in degraded habitats. 


5 Conclusions 


In this study, we used the MaxEnt as a tool for determining the priorities of protecting certain 
endemic species of the genus Onosma in Iran. It seems that the use of this tool can be able to 
accurately predict the potential habitats of these taxa. According to the results of this study, the 
habitat area of these species will increase in the future. Therefore, it is essential to manage and 
conserve the projected areas for these species in the future. 

For successful conservation of plants facing climate change, we recommend regular monitoring 
of all ecosystems in Iran. Plants can be used as indicators at environmental, community, 
population, and individual levels. Thus, we encourage regular monitoring of plants to detect how 
ecosystems are responding to changes. Observation can be done by investigating the presence, 
absence, richness, and composition of plants in defined monitoring stations. Through regular 
monitoring programs, the presence of new species or the absence of previously recorded species 
will be detectable. Moreover, it would be possible to track differences in species distribution, 
richness, and composition. Having this information can support conservation planning of plants 
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dealing with climate change. Regardless of the modeling results, the long-term physiological 
dormancy of seeds and the small populations of some local endemic species, along with the above 
factors may increase the threat of extinction of these valuable endemic species. In addition, 
frequent wildfires and destruction of habitats have put the Zagrosian ecosystems at risk of being 
completely destroyed. The results of this study could justify ex situ conservation strategies such 
as gene banks, field gene banks, and in vitro conservation for these valuable species in Iran. 
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